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ABSTRACT 


This report documents the accomplishments of the NASA sponsored 
effort to develop performance- eva I uat ion capabilities for the design 
of physical systems. These accomplishments are: 

Development of a theory of limiting performance of 
large systems subject to steady state inputs 

. Application and modification of PERFORM, the computa- 
tional capability for the limiting performance of 
systems with transient inputs 

. Demonstration that use of an inherently smooth control 
force for a limiting performance calculation improves 
the system identification phase of the design process 
for physical systems subjected to transient loading. 
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I NTRODUCT I ON 


The primary goals of this study were to formulate and develop a 
capability for the limiting performance of large steady state systems. 

It was shown early in the study that the previously developed capability 
for transient systems could not be extended to steady state environ- 
ments. It was necessary to begin anew and formulate a different theory 
of limiting performance. This new theory has been applied to several 
simple systems. Most of this report deals with the formulation of the 
limiting performance problem for steady state systems. 


LIMITING PERFORMANCE OF STEADY STATE SYSTEMS 

The concept and importance of limiting performance of mechanical 
systems is described in Ref. I. 

The formulation of the limiting performance study of a mechanical 
system subject to steady state sinusoidal loading is given in 
Appendix IA. The formulation is in terms of a linear programming 
problem. The underlying concept in this formulation is the use of 
a Fourier expansion to represent both the responses and control 
forces. The coefficients in the Fourier expansion of control forces 
are the variables to be found. This formulation is applied to a 
single degree of freedom (SDF) system for both linear and nonlinear 
control forces in Appendix IB. The solution for the linear case 
corresponds to results in Ref. I, while for the general case the 
solution compares well with results obtained by an optimal control 
approach (Ref. 2). 

The formulation of Appendix I A is time-dependent. For a linear 
system the problem is reformulated in Appendix II A as a time- independent 
nonlinear programming problem. The merit of this formulation is its 
time-independence with a concomitant reduction in the dimension of the 
problem. The advantage of obtaining the limiting performance of a 
system without multiple analyses of the system has been retained. 
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This formulation is applied to two special cases in Appendix I IB. 

In both cases the results agree with known solutions. 

The steady state limiting performance solution Is extended to 

systems with multiple forcing functions of different frequencies In 
Appendix III. 

The major remaining unresolved problem for the steady state 

formulation is that of limiting performance of systems subject to 

f ore mg functions with frequencies varying over a prescribed range. 

e present formulations are restricted to systems subjected to 

ceding at a prescribed frequency. The study of the frequency range 
problem is continuing. 


APPLICATION AND MODIFICATIONS OF PERFORM 

During the first year of this effort a computer capability 
FORM) was developed for calculating the limiting performance of 
systems with transient inputs. The final report covering this 
effort contains detailed documentation and applications of PERFORM. 

This report has been revised and Is now entitled PERFORM - A PERFORMANCE 
OPTIMIZING COMPUTER PROGRAM FOR DYNAMIC SYSTEMS SUBJECT TO TRANSIENT 
LOADINGS. The revised report contains new results for a train impaot 
problem and a STOL ride control problem. The application to the train 
impact problem was used as part of the paper "Limiting Performance of 
Ground Transportation Vehicles Subject to Transient Loading" (Ref ,, 
presented at the AIAA/ASME/ASCE 13th Structures, Structural Dynamics 

Ref 7;:rrr- ,n +he probiem - the tra,n ,mpac+ 

Ref. ,s modified so that the cushion or shock absorber Is replaced 
by a control force. The problem is to find the minimum force that 
must be transmitted to the lading for a given cushion travel distance 
under specified Impact conditions. Performance tradeoff diagrams were 
rnpu e . tails of the problem formulation are given in Ref. 3 . 
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The limiting performance problem of the ride control system 
for the Twin-Otter STOL airplane (Ref. 5) has been put in PERFORM 
format. Computations of the sort* shewn- i-nr fig'. I were made. Here 
the tradeoff between the min-max acceleration at the center of gravity 
of the airplane and the level of controls was calculated using 
PERFORM. Details of the formulation are to be found in Ref. 5. 

PERFORM has been modified to permit its use in conjunction with 
OPTIMA, a CDC computer software system that can solve large linear 
programing problems. Previously, the IBM program MPS/360 was used. 

SYSTEM IDENTIFICATION 

The use of a Fourier series expansion of control forces in the 
steady state problem suggested a possible application for the design 
of transient systems. It was shown that use of the "smooth" control 
force employed in the steady state problem for the limiting performance 
of transient systems eases the task of Identifying the corresponding 
suboptima I design configuration. Heretofore, transient systems had 
been designed using the limiting performance based on control forces 
that were given a piecewise constant time discretization. This 
application of the Fourier series to limiting performance problems 
was demonstrated for single and two degree of freedom systems. 
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APPENDIX I A 


Formulation of the Limiting Perfor mance Proble m for a Dynamic 
System Subject to Steady-State Sinusoidal Disturbances . 

Consider a dynamic system described by the following equations 
of motion: 


MX + 

• 

CX 

+ 

'3. 

+ 

ii 

Ff 

(1 ) 

where M_ 

= 

N x N 

mass matrix 


C 

= 

N x N 

damping matrix 


K 

= 

N x N 

stiffness matrix 


V 


N x J 

coefficient matrix associated 

with 




control force vector 


F 


N x L 

coefficient matrix associated 

with 




forcing function vector 


and N 

— 

no. of degrees 

of freedom of the system 


J 

= 

no. of isolators or controllers 


L 

= 

no. of disturbances applied on the system 




where 


Assume 

U 


|f L 


m 

E 

i = 


U s jSini(u+ + U c . cos i tot 


(3) 


where 





U, - 
lei 

0 , = 
SI 

U 2si 

• 

Cl 

o 

(1 

o • • 

CNI 

ZD 


• 

“oil 


• 

U,* , 

_ Jc t_ 


(3a) 


m = number of pairs of terms used in the Fourier representation 
of 0 


Simi lari y, assume 
m 

X = E 

i = I 


X . s i n i tot + 
s i 



cositot 


where 



" X lsi" 


o 

_Z i 

X . = 

SI 

X 2s i 

• 

• 

• 

XI 

o 

II 

X 2ci 

• 

• 

• 


. X NsL 


• 

X K, t 

Nc L) 


(4) 
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Def i ne 



x . 

S 1 


0 . 
SI 

X 

II 

X . 
ci _ 

0. = 
1 

1 

Cl 1 
O 1 

— 1 

1 


Substitute (3), (4) into (I) and use (5) to obtain 


X. = 


W.U. + Z.f_ 

—ii —i 0 


where 


W. 

—i 


Y. 

—i 


-1 


i 

-V i 

i 

1— 

l 

0 I 

— I 


0 


wss. ! 

wsc:. 



i 1 
1 - 

1 

-V 


1 

WCS. I 
1 1 

wcc. 



■« m f“ 


“zs ' 

= Y .- 1 

1 1 — 

II 

—i 

0 


ZCj 


_ _ 



6 li = 


I if i = 
0 if i ^ 


Y. = 

—i 


-Mu 2 i 2 + K l 

~ “ I 

h- 


-Ciu 


CIoj 


j -Moj 2 i 2 + K 


( 5 ) 


(6) 


(7) 


Thus, from equation (6), all the coefficients in the Fourier ex- 
pansion of the responses (i.e, X’s) can be expressed in terms of those 
of U's. 


From Eqs. (6) and (4) we get 


X 


= R U + T 
— m m — m 



( 8 ) 
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where 


R 

-m 


WSS_j si not + WCSj cosot WSC ^ si not + WCC ^ cosot 

•WSS sinmot + V/CS cosmot WSC sinmot + WCC cosmot 

m m m m 


m 

T = I 
" m |~| 


ZS . siniot + ZCj cosiot 


(9) 


U = 
m 


m 


and from (3) 

0 = SC 0 


where 


si 

0 . 

C I 


sm 


cm 


SC = 

— m 

I = J x J identity matrix 


I si not I co suit | 
I “ I 


I i 

i I sinmot i I cosmot 

I _ i — 


(10) 


We have expressed X, 0 in terms of 0 , a set of unknown 

m 

numbers. Note that the matrices T , R , SC are all functions of time. 

— m — fn — m 

Suppose the limiting performance problem of interest can be stated 
as: Find the 0 that minimizes the maximum ? while subject to the con- 
straints 


YL. * 0, * YU. 
i i i 


where 



( 12 ) 


__ «* • 

0 = Yj_ X + Y2_ X + Y3 X + Y4_ 0 + Y5^ f 

This problem can be converted to: Find <f>, 0 such that <j> 

subject to 

~ <t> for j = I , 2,. . .NOB 

YL. < 0. ^ YU. for i = I, 2,...N0C 

where 

NOB = number of objective functions 
NOC = number of constraints 

The problem defined by (13) is one of linear programming, i.e. 
Find z 

to minimize c z 

> 

subject to H z = G 
~ < 

where 

£ =[l, 0, 0..., o] 


-T 

> RG 

_ 


+T 

1 RG 
i — m 



0 

! RY 


i — m 


— i 

0 

! RY 

• 

! — m 



where 


G = 


-TG 

f A 

— m 

0 

-TG 

f A 

— m 

_ 0 

YU - 

TY K 


— m 0 

YL - 

TY f A 

. 

— m 0 


RG 
— m 


fi + P2 Rl + P3 R + P4 SC 
m m rn 


i s minimized 


(13) 
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TG 

= PIT2 

+ 

P2TI 

+ 

P3T 

+ P5 f si nwt 

— m 

m 


m 


— -m 

— o 

RY 

= YIR2 

+ 

Y2RI 

+ 

tr 

to 

>- 

+ Y4SC 

— m 

m 


m 


m 

ffl’ 

TY 

= YIT2 

+ 

Y2TI 

+ 

Y3T 

+ Y5f si nut 

— m 

m 


— —fn 


— -m 

— o 


Here R , T , are functions of time and 
— m — m 



T2 

— m 


Tl 


— m 


= R (t) 

-tn 

= R (t) 

-m 

• • 

= T (t) 

-m 

= T (t) 

-m 


li 



APPENDIX IB 


\ 


Solutions of Limiting Performance Problem for Single Degree of Freedom 
System 


Consider the SDF problem (Fig. 1-1) of finding U such that 
cp = max |Z| = max |U| is minimized while jxj - A 


The equation of motion Is 
mZ + U = 0 

Z = X+ f = X + sl nwt 
These give 

X t U = w 2 sinoit 
for a unit mass. 



For A = 0.5, w = 20, the optimal control approach (Ref. 2) gives mlmj> = 162.4. 
The same problem has been solved by linear programs! ng using 9 pairs of 
terms in a Fourier representa lion of U, and the result Is min<J> = 163.2. 

These solutions are very close. The optimal control approach is prac- 
tical for simple systems only, while the linear programming formulation 
applies to large systems. 

For the first time in limiting performance studies it is possible 
to distinguish between linear and nonlinear control forces. If the 
Fourier series Is restricted to one sine and one cosine term then the per- 
formance of a linear system is obtained. For the linear SDF problem, 
the linear programming formulation provides the same solution calculated 
for this case in Ref. I. The new formulation Is significant because the 
treatment of Ref. I cannot be extended to complex systems to which the 
linear programming formulation, can, be app I ied. 
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APPENDIX I I A 


Formulation of a Time-Independent Limiting Performance Solution for Linear 
Dynamic System s Subject to Sinusoidal Inputs 

Consider a multidegree of freedom system, described by 

MX t CX + KX + VU = £_ f o s I mot ( 1 5 


Since we seek the limiting performance for a linear system, let 

0=0 sinwt + 0 coswt 
s c 

X = X si nwt + X coswt 
s c 

Substitute (2) into (I) to give 


V 

= w 

0 

s 

s. 


1 

o 


( 2 ) 


(3) 


where _l 


r — 


r i 


_ — 

^w 2 + K -Cw 


ol 

: T I 


wss wsc 

Cw -Mu 2 + K 


0 -V 


wcs wcc 

__ — 


- - 


— - 


-1 


z = 

-Mw 2 + K -Cw 


F_ 


■ " 
zs 





. 



Cw — Mw 2 + K 


0 


ZC 

— — 


Let the objective function be 

? = PIX t P2X + P3X + P4U + P5f (4) 


and the constraint be 



Using (2) and (3), (4) and (5) beeofZ© 


p = $ 


s 


s i nut + <f 

c 


cos tut 


and 

YL * 0 si ncot + 0 cosut; ** YU 
s c 

where 

= -Pl_ WSS_ U s cu 2 -PJ_ WSC_ U c u 2 -PJ_ ZS_ ? 0 o) 2 
-P2 WCS 0 s u> -P2 WCC U c w + P3 WSS_ U s 


(6) 

(7) 

(8) 


+P3 U c + — — ^0 + — + — ? o 

The quantity $ c is the same as $ s with WCS , V.'CC replacing WSS , WSC re- 
spectively. 0 and 0 are the same 1 as £ , £ with PI P5 replaced 

r 'sc sc 

by Yl Y5. 

From (6), we can write 

max Y. = $ . 2 + <J> , 2 (9) 

i si cl 

and (7) is equivalent to 

0. 2 + 0. 2i YU. 2 (10) 

si ci i 

since, in general, YL ? = -YU. 

Hence, if we define 

d> = max O . 2 + $ . 2 1 for I = 1 , 2 , < • • , NOB 

i 

the limiting performance problem can be restated as: Find U g , U c 

such that 


<j> is mi n imi zed and 


$ 2 t $ . 2 - (J> for i = I, 2, NOB 

si ci 

and 


(II ) 


0 . 2 + 0 . 2 < YU 2 j = I, 2, NOC 

sj cj j 
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CZI 


where NOB = no. of objective functions 
NOC = no. of constraints 

Using (8), everything in (II) can ; he-e^iressfed? i n terms of 0 g and 
. Note that the time dependency of the response variables has been 
eliminated. The problem posed by (II) is one of nonlinear programming. 
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APPENDIX I IB 



Formulation 


|. SDF system (Fig. I 1-1) 

The equations of motion are 
mZ + U = 0 

Z = X + f = X + f Q slnwt 

If m = I 
• • 

X + U = u> 2 f 0 slnwt 

Figure I I -I ( I ) 

We want to find U that minimizes |u| subject to |x| -A 
Let 

U = U s sinwt + U c cosuit (2) 

X = X sinwt + X coswt 
s c 

Then, following the previous formulation. 




The problem now is to find U that minimizes (U,. 2 + U c 2 ) subject to 


-S + ( _£.)2 - A 2 


(3) 


The solution can be found easily (e.g. by graphical means) 
We get 

U = w 2 (f n - A) 
s U 


max|u| = w 2 (fg ~ A) 
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Equation (4), when normalized, is the same as the one obtained In Ref. I 
2. Two Degree of Freedom (DOF) System 

For the two DOF system shown in Fig. 1 1 —2a , Den Hartog (Ref. 6) 
derived the result that under certain 
conditions, the main mass M In the 
system does not move at all. 

This optimal condition can be obtained 
using the limiting performance approach. 

The equations of motion are (Fig. Il-2b) 

MX x + KXj + U = P Q s I nut 



mX2 - U = 0 


Then 


U = U si nut + U coswt 
s c 


X I = X si slnu)+ + X ci cosw+ 1 = 1 * 2 

ul-Pn “ U 

X. 

Mu 2 -K *■'' mu' 

-U 


I s 


_s 

.,2 



U 


Me Mu 2 -K 


X 2c~ 


mu 


(b) 

Figure 11-2 


Now, we want to find the U that minimizes | Xx J . This Is equivalent to 


finding the U that minimizes 


X 2 + X 2 = 
is lc 


U -P n 
s 0 


U 


Mu 2 -K 


^Mu 2 -K 

The solution is seen to be 

U = p„, (J = 0 and min maxlXjl = 0 
s 0 c 

Suppose now we modified the problem by putting a constraint on the rattle 

space between the two masses. Now we seek min max|Xx| while |Xx-X 2 | 5 A. 

Stated In terms of a nonlinear programming problem, this becomes: Find 

U , U to minimize 
s' c 
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X. 2 
Is 


+ X 


lc 


subject to 

<X 1S - X 2 s )2 + <X lc- X 2c )2 A2 
The solution (valid for the case Mu 2 > K) is 


U s = P 1 + T,,, U c = ° 


min max|Xi = 


P ! tx l - P C 
|Mw 2 - k| 


where 


P. = 


P 0 mw 2 


■ -Mu 2 + K + mus 2 


Y = 


mw 2 (Mo) 2 -K) 


Mu 2 -K-mw 2 


This problem would be very difficult If it were to be treated using the 
Den Hartog approach. 
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where 


PRECEDING PAGE BLANK NOT FILMED 



Substitute (2), (3) into (I) and equate the coefficients of terms 

coslw.t and siniu t respectively. Then 
J J 


(-Mi 2 uK 2 + K) X js| (Cl a) j ) Xj cj Fj 4 (1 — ^jsi 

(CiMj) X jsi (-Mi 2^. 2 + K) X jc! = lO jci 



or 
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and 


Z J ' = 


■J ' 


zs 


zc ji 


= Y 


lj 


,U _ 


-Ml ^u . ^ + K 
- J 


L -“j 


0 


-C| U j 


-T -1 


-M i ^u . ^ + K 
- J 


From (3) and (4) we have 


L m 
X = Z I 

j=l i=l 


R j 1 0 j ' + f j 1 

-mm m 


(5) 


where 


R 

-m 


J ' _ 


WSS J ‘ * siniu.t + WCS J 1 coslu.t 
J J 


WSC J * siniu.t + WCC J ' cosiu.t 


T Jl = 
m 


ZS^ siniu.t + ZC^ ' coslu.t 
J J 


a J' = 

m 


U jsi 

lVu 


Equation (2) can be rewritten as 


L m j ; = : | 

U = Z Z SC J U J 

j = | | = | “* m m 


(6) 


where 


SC 

— m 


ji _ 


[ 


| siniu.t \ | cosiu 
— J • 


>•] 


I is a J x J identity matrix. 
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We have just expressed X and 0 in terms of the unknown coefficients 
0 j». Following the procedure In Appendix IA, the limiting performance 
problem can now be placed in linear programming form with the coefficients 

of 0 as unknowns. 
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